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ABSTRACT 



Context. Models have been made of stars of a given mass which produce planetary nebulae which usually begin on the AGB (although 
they may begin earlier) and run to the white dwarf stage. While these models cover the so-called dredge-up phases when nuclear 
reactions occur and the newly formed products are brought to the surface, it is important to compare the abundances predicted by the 
models with the abundances actually observed in PNe. 
c I ^ , Aims. The aim of the paper is to determine the abundances in a group of PNe with uniform morphological and kinematic properties. 

04 . The PNe we discuss are circular with rather low-temperature central stars and are rather far from the galactic plane. We discuss the 

effect these abundances have on determining the evolution of the central stars of these PNe. 

Methods. The mid-infrared spectra of the planetary nebulae NGC 1535, NGC6629, He2-108, and Tel (IC 1266) taken with the 
Spitzer Space Telescope are presented. These spectra were combined with the ultraviolet lUE spectra and with the spectra in the 
^0 ' visual wavelength region to obtain complete, extinction-corrected spectra. The chemical composition of these nebulae is then found 

by directly calculating and adding individual ion abundances. For two of these PNe, we attempted to reproduce the observed spectrum 
by making a model nebula. This proved impossible for one of the nebulae and the reason for this is discussed. The resulting abundances 
are more accurate than earlier studies for several reasons, the most important is that inclusion of the far infrared spectra increases the 

Q . number of observed ions and makes it possible to include the nebular temperature gradient in the abundance calculations. 

Jh ' Results. The abundances of the above four PNe have been determined and compared to the abundances found in five other PNe with 

tZt. ' similar properties studied earlier These abundances are further compared with values predicted by the models of Karakas (2003). 

j^ , From this comparison we conclude that the central stars of these PNe originally had a low mass, probably between IM© and 2.5Mo. 

A further comparison is made with the stellar evolution models on the HR diagram, from which we conclude that the core mass of 
these PNe is between 0.56Mo and 0.63Mo. 

Conclusions. A consistent picture of the evolution of this group of PNe is found that agrees with the predictions of the models 
concerning the present nebular abundances, the individual masses, and luminosities of these PNe. The distance of these PNe can be 
determined as well. 

(3 . Key words. ISM: abundances - planetary nebulae: individual: NGC 1535 & Tc 1 (& IC 1266) & NGC 6629 & He2-108 - Infrared: 

^H ■ ISM: lines and bands 

,__! 1. Introduction argon, sulfur, chlorine, and sometimes other elements are de- 

__l termined as well. These elements are ejected by the star and are 

; The evolution of planetary nebulae (PNe) has been studied for ^^easured in the nebula. These abundances provide complemen- 

> . about four decades with the result that a general picture of the information which aid in understanding the evolution and 

r^ evolution is understood, but the details are still being debated. ^^^-^^^ ^^„ eventually confirm or limit the results of models of 

rN This is both because many of the parameters involved are poorly j^ evolution 

S known and because the atmosphere of the central star is difficult 

■ ■ ■ to describe, including the eff^ective temperature, radius, and mass We are concerned here with abundances found in PNe whose 

loss rate. The (in many cases) uncertain distance enhances these central stars are rather bright and whose morphology is usually 

difficulties. Furthermore, the interaction of the central star with described as round or elliptical. We present new abundances for 

the nebula is often unclear. This manifests itself in a difficulty ^ur such nebulae: NGC1535, NGC6629, Tel (IC 1266), and 

explaining the ionization structure in some PNe. He2-108. The abundances in these nebulae will be combined 

The abundances of some elements change in the course of with the results found earlier in five other PNe, which ai-e also 

evolution of the central star. We discuss here primarily the el- excited by bright stars and which have similar (round or elhp- 

ements helium, nitrogen, carbon, and oxygen, although neon, ^ical) morphology. These PNe have similar spatial properties as 

well: they are at rather high galactic latitudes. We can therefore 

Send ojfprint requests to: pottasch@astro.rug.nl speak of a class of nebulae. This class has been extensively stud- 

* Based on observations with the Spitzer Space Telescope, which ied earlier but from a different point of view. The central stai's 

is operated by the Jet Propulsion Laboratory, California Institute of of these PNe are among the brightest known in the visible. It is 

Technology. therefore possible to observe these spectra with very high res- 

** formerly with the Indian Institute of Astrophysics, Bangalore olution. The profiles of hydrogen and helium lines can then be 
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analyzed with the initial goal of obtaining the effective temper- 
ature (Teff) and gravity of these stars. The final goal is to use 
these quantities to study the mass and luminosity (thus the evo- 
lution) of these stars, but additional information is required to do 
this. 

One of the early ex t ensive studies of the line profiles is that 
of 'Mend ez et all (Il988l Il992h . These authors take high resolu- 
tion optical spectra of 24 PNe and interpret these spectra us- 
ing NLTE plane-parallel static model atmospheres which con- 
tain only hydrogen and helium. The values of T^/f and the stel- 
lar helium abundance are determined by fitting the observed 
Hel and Hell line profiles while the gravity (log g) is found 
by fitting the hydrogen line profiles (usually Hy but some- 
times HyS). Then, making use of theoretical evolution diagrams 
which plot Tgff against log g for different values of the stellar 
mass M5/M0, the value of M.5/M0 is determined. Using this 
mass and the value of gravity, the angular radius of the ex- 
citing star is determined. Combining this value with the pre- 
viously found temperature and the observed magnitude cor- 
rected for extinction, the distance to the nebulae is found. In 
the course of time this group has improved the model atmo- 
spheres used by first i ncluding the sphericity of the atmosphere 
dKudritzki et al.|[l997h and then by including ot her elements be- 
sides H and He in the atmosphere dKudritzki et al. 2006). These 
improved models do not change the stellar masses or PNe dis- 
tances substantially. Both the masses and the distances deter- 
mined in this way, however, are suspect. Quoting Kudritz ki et alJ 
d2006l) 'the masses determined seem systematically larger than 
white dwarf masses and some of the objects ... have unrealisti- 
cally high masses'. Concerning the distances: for the 18 cases 
where the di stance thus found may be co mpared with statisti- 
cal distance dStanghelUni et alj |2008: Cah net alj|1992h they are 
always highe r, in 11 cases more than 60% higher 

Recently iPauldrach et alj d2004l) used a different approach. 
By including models losing mass they are able to use 'ob- 
served' mass loss rates and terminal velocities to obtain the stel- 
lar masses and distanc es. These masses and distances are even 
higher than found by Kudritzki et al. (2006') and are shown to 
be implausibly high by Napiwotzki (2006). Two arguments are 
used. Because high mass central stars are probably descended 
from high mass progenitors they should have the properties of 
high mass stars. First they should belong to the disk population 
of the galaxy which is characterized by small scale heights per- 
pendicular to the galactic disk. Secondly such high mass star 
should have produced and dredged up substantial amounts of 
both nitrogen and helium. Napiwotzki conclud es that on bot h of 
these points the high masses found by Pauldrac h et al.l d2004l) are 
unlikely. 

This suggests the following step: the abundances of those 
elements produced by the various stars should be accurately 
determined. These abundances may then be used to determine 
the mass of the star in question by comparing the observed 
abundance with the abundances determined by nucleosynthesis 
model calculations made in the course of evolution of stars of 
different masses. For these models, the calculations made by 
iKarakasI d200 3') are used. These masses can then be compared 
with those found by Kudritzki et al. (2006). Because the masses 
used in the calculations o f Karakas are th e initial stellar masses, 
and the masses given by iKudritzki et al.l (12006 ) are the present 
(almost final) stellar masses, a ratio between the initial and final 
mass must be known. For this secondary information the values 
given by Weidemann (2000) will be used. Unfortunately the un- 
certainties present in each of these steps will accumulate so that 
definite conclusions are difficult to draw. It is however important 



that central star masses can be derived from accurate nebular 
abundances. 

The main purpose of this paper is to obtain accurate abun- 
dances for four of the nebulae discussed by ?. The most impor- 
tant reason that this can be achieved is the inclusion of the mid- 
infrared spectrum taken with the IRS spectrograph of the Spitzer 
Space Telescope (Werner et al. 2004). The reasons for th is have 
been discussed in earlier papers (e.g. see Pottasch & Be intemal 
ll999l:lPottaschetal . 2000, 2001; Bernard-Salas et al. 2000), and 
can be summarized as follows: 1) The intensity of the infrared 
lines is not very sensitive to the electron temperature nor to pos- 
sible extinction effects. 2) Use of the infrared line intensities en- 
able a more accurate determination of the electron temperature 
for use with the visual and ultraviolet lines. 3) The number of 
observed ionization stages is doubled. 

For two of the PNe a second method of determining the 
abundances is attempted using a nebular model. This has sev- 
eral advantages. First it provides a physical basis for the electron 
temperature determination. Secondly it permits abundance de- 
termination for elements which are observed in only one, or a 
limited number of ionic stages. This is true of Mg, Fe, P and CI, 
which would be less reliably determined without a model. A fur- 
ther advantage of modeling is that it provides information on the 
central star and other properties of the nebula. 

A disadvantage of modeling is that there are more unknowns 
than observations and some assumptions must be made espe- 
cially concerning the geometry and the form of the radiation field 
of the exciting star. The difficulties we encountered in determin- 
ing a model for NGC 1535 are discussed below. 

Abundance determination for these nebulae have been made 
earlier and will be compared with our values in the individual 
sections. 

This paper is structured as follows. In Sect. 2, 3, 4 and 5 the 
observations of NGC 1535, NGC 6629, Tel and He2- 108 will be 
discussed and the resultant abundances will be given and com- 
pared to earlier results. In Sect. 7, the models for NGC 1535 and 
Tel are presented and discussed. In Sect. 8 the evolutionary state 
of the nebulae is discussed. Finally, our conclusions are given in 
Sect. 9. 



2. NGC 1535 

NGC 1535 (PN G206.4-40.5) is a roughly circular nebula con- 
sisting of a bright inner ring of about 20 "x 17"surrounded b y 
a much fainter outer shell of 48"x 42" (.Banieriee et al]|1991h . 
iTvlenda et all d2003h list a size of 33"x 32"down to the 10% 
level. The nebula is at a very high galactic latitude. The distance 
to the nebula is rather uncertain. Ciardullo et al. (1999) give a 
value of 2.3 kpc o n the basis of a p ossible association of a nearby 
star with the PN. iHerald & Bianchi (.2004.) use a value of 1.6 
kpc on the basis of nebular models. We shall use the latter value 
when necessary for two reasons. First we feel that the evidence 
for association of the star and nebula is rather weak. Secondly at 
the larger distance the dimensions of the nebula are rather large 
for such a bright nebula. We stress however that the distance is 
uncertain. 



2. 1 . The infrared spectrum 

Observations of N GC 1535 we re made using the Infrared 
Spectrograph (IRS iHouck et al.l 2004) on board the Spitzer 
Space Telescope with AOR keys of 4111616 (on target) and 
4111872 (background). The reduction started from the droop 
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images which are equivalent to the most commonly used Basic 
Calibrated Data (bed) images but lack stray-cross removal and 
flat-field. The data were processed using th e si 5. 3 version of th e 
pipeline and using a script version of Smart dHigdon et alJ20041) . 
The tool irsclean was used to remove rogue pixels. The different 
cycles for a given module were combined to enhance the S/N. 
Then the resulting high-resolution modules (HR) were extracted 
using full aperture measurements, and the SL measurements set- 
ting a window column extraction. This same reduction has been 
used for all nebulae discussed in this paper. 

The IRS high resolution spectra have a spectral resolution of 
about 600, which is a factor of between 2 and 5 less than the 
resolution of the ISO SWS spectra. The mid-infrared measure- 
ments are made with several diff'erent diaphragm sizes. Because 
the diaphragms are smaller than the size of the nebulae and are 
all of differing size, we first discuss how the different spectra are 
placed on a common scale. 

Two of the three diaphragms used have high resolution: the 
short high module (SH) measures from 9.9^(m to 19.6;um and the 
long high module (LH) from 18.7/im to 37.2yum. The SH has a 
diaphragm size of 4.7"x 11.3", while the LH is ll.l"x 22.3". 
If the nebulae are uniformly illuminating then the ratio of the 
intensities would simply be the ratio of the areas measured by 
the two diaphragms. Since this is not so, we may use the ratio 
of the continuum intensity in the region of wavelength overlap 
at 19;um. These continua are equal when the SH intensities are 
increased by a factor of 3.2. The third diaphragm is a long slit 
which is 4" wide and extends over the entire nebula. This SL 
module measures in low resolution and measures between 5. 5;um 
and 14yL(m. These spectra are normalized by making the lines in 
common between the SL and SH modules to agree. Especially 
important is the agreement of the [S iv] line at 10.5 lyum. A uni- 
form scale is obtained by increasing the SH intensities by a factor 
of 3.2 with respect to the LH intensities. The SL intensities are 
increased by a factor of 1.23. 

The IRS measurement of NGC 1535 was centered at 
RA(2000) 04''14'"15.9" and Dec(2000) -12°12'21 ". This is al- 
most e xactly the same as the value measured by i Kerber et alj 
dlOOl of RA(2000) 04''14"'15.78' andDec(2000) -12°44'21.7", 
which is presumably the coordinate of the central star Thus the 
IRS measurement was well centered on the nebula. The fluxes 
were measured using the Gaussian line-fitting routine. The mea- 
sured emission line intensities are given in Table 1, after cor- 
recting the SH measurements by the factor 3.2 and the SL mea- 
surements by a factor of 1.23, in the column labeled 'inten- 
sity'. The Hj6 flux found from the infrared hydrogen lines (es- 
pecially the lines at 12.3 7^^m) using the theoretical ratios of 
iHummer & StorevI (Il987h . is 2.3x10 "erg cm"^ s ', which is 
about 49% of the total Hj6 intensity. This is reasonable since 
the well-centered LH diaphragm covers a large fraction of the 
nebula. Note that by scaling in this way the nebula is assumed 
homogeneous when it is larger than the different modules. 

2.2. The visual spectrum 

The visual spectrum has been measured by at least six authors. 
We list here the results from four of these. The line intensities 
listed have been corrected by each author for a value of extinc- 
tion determined by them to obtain a theoretically correct B aimer 
decrement. The result are listed in Table 2, where the last column 
lists the average value which we have used. No attempt has been 
made to use a common extinction correction because then the 
Balmer decrement will be incorrect. The value of extinction C 
which the individual authors found is listed at the bottom of the 



table. None of the spectra measure the weaker Unes very well. 
The errors may be judged by the agreement (or disagreement) 
of the various measures and appear to be within 20% for the 
stronger lines and worse for the weaker lines. 



Table 2. Visual spectrum of NGC 1535. 



Ion 



Intensities 



Average 



(A) 




(1) 


(2) 


(3) 


(4) 


In tens. 


3727* 


[On] 


4.1 


7.1 


8.41 


9.37 


8.41 


3869 


[Nem] 


90 


95.8 


97 


116 


101 


4267 


Cii 


0.4: 


0.34 


1.03 


0.36 


0.34 


4363 


[Ora] 


13.3 


12.5 


11 


12.6 


12.6 


4686 


Hen 


27 


14.2 


18.3 


17 


17.5 


4711 


[Any] 


4.5 


4.2 


5.4 


4.5 


4.5 


4740 


[Any] 


3.8 


3.2 


3.8 


3.36 


3.4 


4861 


W 


100 


100 


100 


100 


100 


5007 


[Om] 


1220 


1210 


1180 


1190 


1200 


5517 


[Clm] 




0.242 


0.264 




0.25 


5538 


[Clm] 




0.177 


0.175 




0.175 


5755 


[Nn] 




0.1: 


0.62 




0.1: 


5876 


He I 


10.2 


12.6 


11.8 


12.6 


12.2 


6312 


[Sm] 




0.303 




0.28 


0.30 


6584 


[Nn] 


1.8 


1.62 


0.25: 


0.85 


1.6: 


6717 


[Sn] 




0.0627 






0.0627 


6731 


[Sn] 




0.114 






0.114 


7135 


[Arm] 


6.0 


6.42 


5.3 




6.3 


7263 


[Any] 




0.217 






0.217 


8045 


[CI IV] 




0.512 


0.53 




0.52 


9532 


[Sm] 




6.75 






6.75 


C(H0) 




0.2 


0.07 


0.11 


0.01 





t References; (1) Barker ( 1989), (2' Milingo et"ai] ( l2010h . (3) lAller et afl 

( 1979), (4) ,Ki-abbe & CopettL ^2006) . 

(:) indicates uncertain values. 

* This is a blend of A3126 and A3129 lines. Only Aller & Czyzak are 

able to resolve this doublet: 3726=5.4, 3729=3.01 

C is the extinction used by the author 



2.3. The ultraviolet spectrum 

Quite a large number of lUE spectra of this nebula have 
been taken. If only those taken with the large diaphragm are 
counted there are six shortwavelength spectra and five long- 
wavelength spectra with the central star included. These have 
been made with low resolution. Also two shortwavelength spec- 
tra (SWP 10821 and 13495) and one longwavelength spectrum 
(LWR02165) were taken with high resolution. They also in- 
cluded the central star Two spectra are also available which 
included only the nebula (SWP15497 and LWR11975). These 
spectra are especially useful because longer exposures can be 
made without saturating the stronger lines. The N in line can only 
be seen on these spectra. The lUE diaphragm is an ellipse about 
10"x 21"; the wavelength range is from 1 150A to about 3220A. 
The high resolution spectra have a resolution of 0.2A while the 
low resolution is about 6A. There are also two spectra taken 
with the Hopkins Ultraviolet Telescope (HUT), one taken with 
the central star in the diaphragm, the other has only the nebula 
in the 9.4"x 116"diaphragm. The HUT spectra have a resolu- 
tion of about 3A and they cover the wavelength range between 
830Aandl860A. 

Because the nebula is in all cases larger than the diaphragm 
used and the various spectra are centered at different positions 
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Table 1. IRS spectra of the four nebulae. The measured line intensity is given in Col. 3, 5, 7 and 9. Col. 4, 6, 8 and 10 give the ratio 
of the Hne intensity to HyS(=100). 







NGC 1535 


NGC 1535 


Tel 


Tel 


He2-108 


He2-108 


NGC6629 


NGC6629 


Identification 


A(fim) 


Intensity^ 


1/^8 


Intensity" 


I/H/3 


Intensity'' 


I/H/? 


Intensity" 


I/HyS 


F 


6.48 






198±27 


3.2 










[Nin]? 


6.63 






76.5±8 


1.24 










[Am] 


7.026 






2020±70 


32.8 










Hi (6-5) 


7.47 


36.2+2.1 




230±17 




40.7±4.8 




306±29 




F 


7.757 














37.2±13 


0.44 


F 


8.497 






384±10 












[Arm] 


8.99 


76.5±4.3 


3.32 


399±16 


6.50 


299±13 


25 


1080±45 




[Siv] 


10.512 


915±35 


39.8 


29.2±2.2 


0.47 


41.0±2.4 


3.6 


1220±46 


14.4 


Hi (9-7) 


11.305 






20.6±3 








32.4±2.3 




[CI IV] 


11.763 


10.5±1.9 


0.455 














Hi (7-6-1-11-8) 


12.375 


22.4:±2.1 




67.9±3.6 




11.7±1.2 




95.8±4.3 




[Neil] 


12.81 


9.55±1.9 


0.415 


2301±200 


37.5 


1175±130 


102 


839±19 


102 


[Neiii] 


15.555 


1920±40 


83.5 


89.5±5 


1.46 


189±6 


16.5 


7530± 


89.1 


F 


17.66 










35±12 


3 






[Piii] 


17.88 














51±2.3 


0.61 


[Siii] 


18.714 


67.4±3.3 


2.93 


863±80 


14.0 


830±150 


72.2 


1540±32 


18.2 


[CI IV] 


20.316 


10.8 ±0.9 


0.47 














[Arm] 


21.820 


7.44 ±0.73 


0.323 


26.6±2 


0.432 


19.4±2 


1.69 


85.4±6.9 


1.01 


[Fem] 


22.92 






18.1±1.6 


0.295 


11.7±0.7 


1.02 






[Oiv] 


25.889 


2466 ±21 


107 














[Siii] 


33.48 


63.2 ±1.9 


2.75 


382±33 


6.21 


402±18 


35.0 


807±34 


6.3 


[Nem] 


36.013 


203±2.2 


8.84 










531±24 


6.3 



t Intensities measured in units of 10 '''erg cm ^ s ' 



The intensities SH measurements (below 19yum) and the SL measurements have been 
increased by a factor to b ring them all to th e scale of the LH measurements. The factors used are given in text. Lines identified by T' are related 
to the FuUerine molecule dCami et al.ll2010h . 



in the nebula, a total spectrum is obtained by normalizing the 
individual spectra to the strong Hell line at 1640A. Thus as 
basis we use the high resolution lUE spectra SWP 13495 and 
LWR02165; the other lines are found by using the measured 
ratio to the intensity of the 1640A line in the individual spec- 
tra. The values are given in the third column of Table 3. The 
measured values are then corrected for the diaphragm size and 

erg cm" 



Table 3. UV Spectrum of NGC 1535. 



the extinction using the total H/3 flux of 4.73 x 10 " °'"" ""^ ^ 



s" (see below), then the intensity of the Hen at 4686 A be- 
comes 8.25 X 10"" erg cm"^ (from Table 2). The theoretical 
helium spectrum (Hummer & Storev 1987) then gives the inten- 
sity of the 1640 A line to be 5.64 x 10~" erg cm"^. The other 
lines are then corrected for their extinction with respect to the 
1640A line using the values given by Fluks et al. (J,994) and the 
value C=0.09 as given below. These values are listed in col. 4 
of Table 3. In col. 5 of the table the ratio of the line to H/3 nor- 
malized to Hy6=100 is given. The uncertainties in the intensities 
are estimated to be 20% for the stronger lines and 30% for the 
weaker lines. 



A 


Ion 




Intensities 


(A) 


(1) 


(2) 


(1/1^0) 


977 


Cm] 


9.6 


8.4 


18 


1175 


Cm] 


34 


23.5 


50 


1548 


Civ 


30.8 


19.3 


40.7 


1550 


Civ 


16.8 


10.0 


21.2 


1640 


Hen 


95 


56.4 


119 


1661 


Om] 


1.8 


1.1 


2.3 


1663 


Oiii] 


6.9 


4.1 


8.7 


1678 


7 


4.8 


2.8 


6.0 


1750 


Nm] 


8.5 


5.05 


10.7 


1906 


Cm] 


85.2 


50.9 


108 


1909 


Cm] 


58.0 


34.5 


73 


2297 


Cm] 


16 


9.5 


20 


2734 


Hen 


3.5 


1.9 


4.0 


2837 


[Feiv]? 


1.8 


1.0 


2.1 


3048 


Om 


7.0 


3.6 


7.2 


3134 


Om 


28.0 


14.2 


30 


3204 


Hen 


7.2 


3.7 


7.7 



2.4. Extinction 



(l)Measured intensity from in units of 10 " erg cm ^ s '. 
(2)Intensity corrected for diaphragm size and extinction in units of 
10" '2 erg cm"2 s"' . I/H/3 is normalized to UJ3= 100 



The two methods which can be used for obtaining the extinction 
are: (1) comparison of radio emission with H/? flux, and (2) com- 
parison of observed and theoretical Balmer decrement. The four 
values of the extinction correction C(Hy6) which are found in the 
literature are given in Table 2, and are seen to have a rather large 
range. Let us discuss the radio emission and the H/3 flux. 



2.4.1 . The 6 cm radio emission and tine H/3 flux 

The 6 cm flux density has been measured bv lGrifflth et alj d 19941) 
as 168 mJy. This corresponds to an H/i flux of 4.73 x 10"" erg 
cm"^ s"' using the electron temperature and helium abundance 
give n below. The obse rved Hj6 flux of 3.8 x 10"" erg cm"^ s"' 
(see ICahn et al.lfl992h . This leads to an extinction constant of 
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C(Hy8)=0.094 or ^b-v =0.064. This is quite similar to the values 
found from the Balmer decrement and listed in Table 2. This 
value, together with the extinction curve of Fluks et al. (,1994b . 
has been used in correcting the UV fluxes in Table 3. 



3. Chemical composition of NGC 1535 

The method of analysis is the same as used in the papers cited 
in the introduction. First the electron density and temperature 
as function of the ionization potential are determined. Then the 
ionic abundances are determined, using density and temperature 
appropriate for the ion under consideration, together with Eq. ( 1 ). 
Then the element abundances are found for those elements in 
which a sufficient number of ion abundances have been derived. 



3.1. Electron density 

The ions used to determine A^e are listed in the first column of 
Table 4. The ionization potential required to reach that ionization 
stage, and the wavelengths of the lines used, are given in Cols. 
2 and 3 of the table. Note that the wavelength units are A when 
4 ciphers are given and microns when 3 ciphers are shown. The 
observed ratio of the lines is given in the fourth column; the cor- 
responding A^e is given in the fifth column. The temperature used 
is discussed in the following section, but is unimportant since 
these line ratios are essentially determined by the density. 

The electron density appears to be about 1000 cm"-' although 
the two ions with the lowest ionization potential give a somewhat 
higher value. These values are less well determined because the 
ratios are poorly measured. The density is probably not uniform 
as indicated by the structures seen in the central area of the neb- 
ula, which may contribute to this diff'erence. A density of 1000 
cm~^ is used in the abundance determination in Table 6, but none 
of the abundances listed in the table is sensitive to the density in 
the range shown in the table. 



3.2. Electron temperature 

A number of ions have lines originating from energy levels far 
enough apart that their ratio is sensitive to the electron tempera- 
ture. These are listed in Table 5, which is arranged similarly to 
the previous table. While there is a slight scatter in these values 
there is no clear indication of a temperature gradient as function 
of the ionization potential as has been seen in some other neb- 
ulae. An electron temperature of 12000 K will be used with an 
uncertainty of less than 1 000 K. 



3.3. Ionic and element abundances 

The ionic abundances have been determined using the following 
equation: 



Mo. 

Mn 



ion ., ^ul "^H^ 



'H, 



Mo, 



(1) 



where hoJIviii is the measured intensity of the ionic line com- 
pared to Hy6, M is the density of ionized hydrogen, A^x is the 
wavelength of this line, A^^ is the wavelength of Wfi, q-h^ is the 
effective recombination coefficient for H/?, Aui is the Einstein 
spontaneous transition rate for the line, and NJN\on is the ra- 
tio of the population of the level from which the line originates 
to the total population of the ion. This ratio has been determined 



usually using a five level atom. Sometimes a two level atom was 
sufficient. 

The results are given in Table 6, where the first column lists 
the ion concerned, and the second column the line used for the 
abundance determination. The third column gives the intensity 
of the line used relative to Hj6=100. The fourth column shows 
the ionic abundances, and the fifth column gives the Ionization 
Correction Factor (ICE). This has been determined empirically, 
usually by looking at the ionization potential of the missing ion. 
Notice that the ICE is unity for all elements except for Ar, S and 
CI where it is close to unity. The helium abu ndance has been 
derived with the help o f the theoretical work of lBeniaminI (Il999h 
and Porter et all (l2005h . 

Table 6. Ionic concentrations and chemical abundances in 
NGC 1535. Wavelength in Angstrom for all values of A above 
1000, otherwise in/ym. 



Ion 


A 


Int./HyS 


A'ion/A'p 


ICF 


A'el./A'p 


He+ 


5875 


12.2 


0.076 






He++ 


4686 


17.5 


0.015 


1 


0.091 


C++ 


1909 


181 


1.23(-4) 






C+3 


1548 


61.9 


3.54(-5) 


1 


1.6(-4) 


N+ 


6584 


1.6 


2.0(-7) 






N++ 


1750 


10.7 


3.2(-5) 


1 


3.2(-5) 


0+ 


3727 


84 


1.38(-6) 






0++ 


5007 


1200 


2.4(-4) 






0+3 


25,9 


107 


2.03(-5) 


1 


2.7(-4) 


Ne+ 


12.8 


0.42 


5.2(-7) 






Ne++ 


15.5 


83.5 


4.70(-5) 






Ne++ 


3869 


101 


4.95(-5) 


1.1 


5.4(-5) 


S+ 


6731 


0.114 


3.7(-9) 






S++ 


18.7 


2.93 


2.29(-7) 






S+3 


10.5 


39.8 


8.9(-7) 


1.2 


1.3(-6) 


Ar++ 


8.99 


3.32 


2.83(-7) 






Ar++ 


7135 


6.3 


3.57(-7) 






Ar+3 


4740 


3.4 


5.75(-7) 


1.2 


l.l(-6) 


C1++ 


5538 


0.175 


1.5(-8) 






C1+-' 


11.8 


0.455 


2.24(-8) 






Cl+3 


8045 


0.52 


3.36(-8) 


1.2 


6.0(-8) 



Intensities given with respect to 11/3=100. 
(:) indicates uncertain values 



3.4. Comparison with other determinations 

In Table 7 the present abundances are compared to earlier deter- 
minations. The agreement is usually within a factor of 2, except 
for Chlorine which is difficult to measure. 



4. Tc1 (iC1266) 

Tel (PN G345.2-08.8, also known as IC 1266, SaSt2-16 
and IRAS 17418-4604) is morphologically quite similar to 
NGC 1535. It is r oughly circular and has a size at the 10% level 
of 12.9"x 12.2" (iTvlend a et al."2003'). A som ewhat smaller di- 
ameter (9.6") is given by Acker et al. ( Il992h . The size is small 
enough so that most of the radiation can be measured in the lUE 
diaphragm. The nebula is surrounded by a faint halo which is 
also circular and has a diameter of about 53". 

The 6cm continuum radio flux density has been measured 
by Griffith e^al. ( 1994) as 140 mJy. Milne & Aller ( 1982) have 
measured the 2cm radio ffux density as 130 mJy, which corre- 
sponds to a value of 147 mJy at 6cm. We use an average value of 
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Table 4. Observed Electron density indicators in the nebulae. 



Ion 


loniz. 


Lines 


Obs. Ratio 


Ne (cm-^' 


Obs. Ratio 


Ne (cm-^' 


Obs. Ratio 


Ne (cm-^^ 


Obs. Ratio 


Ne (cm-' 




Pot.(eV) 


Used 


NGC1535 


NGC 1535 


Tel 


Tel 


He2-108 


He2-108 


NGC6629 


NGC6629 


[Sii] 


10.4 


6731/6716 


1.8; 


2500: 


1.59 


2900: 






0.74 


1600 


[On] 


13.6 


3726/3729 


1.8: 


3000: 


1.51 


1900: 






1.63 


2400 


[Sm] 


23.3 


33.5/18.7 


0.94 


700 


0.444 


3200 


0.485 


2200 


0.525 


2200 


[Clm] 


23.8 


5538/5518 


0.70: 


600: 


1.06 


3000: 






0.9: 


1400: 


[Cm] 


24.4 


1906/1909 


1.48 


900 


1.24 


4500 










[Ariv] 


40.7 


4740/4711 


0.76: 


1000: 















(:) indicates uncertain values 



Table 5. Observed Electron temperature indicators in the nebulae. 



Ion 


loniz. 


Lines 


Obs. Ratio 


re(K) 


Obs.Ratio 


TeCK) 


Obs.Ratio 


re(K) 


Obs.Ratio 


TeCK) 




Pot.(eV) 


Used 


NGC1535 


NGC1535 


Tel 


Tel 


He2-108 


He2-108 


NGC6629 


NGC6629 


[Nn] 


14.5 


5577/6584 






0.0114 


9 300 










[Sm] 


23.3 


6312/18.7 


0.102 


13 000 


0.0337 


9100 






0.029: 


8 700: 


[Arm] 


27.6 


7135/8.99 


1.9 


14000 


0.888 


9000 


0.60 


8 000 


1.0 


8 500 


[Arm] 


27.6 


7135/5192 






190 


9000 










[Om] 


35.1 


4363/5007 


0.0105 


11800 


0.00446 


9000 


0.0055 


9 500 


0.0041 


8 700 


[Om] 


35.1 


1663/5007 


0.0092 


10700 














[Nem] 


41.0 


3868/15.5 


1.21 


12 200 






0.285 


8 200 


0.42 


8 900 



Table 7. Comparison of abundances in NGC 1535. 



Ele. 



pres 



bark 



milin 



AC 



TPP 



KC 



He 


0.091 


0.097 


0.0.96 


0.094 


C(-4) 


1.6 


0.8 


1.91 


3.65 


N(-5) 


3.2 


4.3 


2.09 




0(-4) 


2.7 


3.3 


2.98 


4.05 


S(-6) 


1.3 








Ar(-6) 


1.1 


1.2 


0.937 


1.93 


Ne(-5) 


5.3 


7.1 


6.37 


8.3 


Cl(-7) 


0.60 




0.184 


1.26 



0.091 0.105 



3.8 



9.5 



1.3 

2.6 



7.3 



References: bark: ' Barkej (Il989h. milin: iMilingo et al.l (120 lOl) AC: 
AUeretal. (1979) , TPP : iTorres-Peimbert & Peimb'err ( ll977h . KC: 
krabbe & Copettii ( i2006h 



145 mJy at 6cm, which corresponds to an li/3 flux of 5. 1 x 10" 
erg cm~'^ s~' . The mea sured H/? flux is 2.18 x 10"" erg cm"^ s"' 
fsee lAcker et al.lfl992h which leads to an extinction coefficient 
C=0.36. 



4.1. Infrared spectrum 

The IRS measurement of Tc 1 was centered at RA(2000) 
17H5'"35.3" and Dec(2000) -46°05'23.3" This is almost ex- 
actly the same as the value measured by iKerber et alJ (l2003h 
of RA(2000) 17''45'"35.3' and Dec(2000) -46°05'23.8", which 
is presumably the coordinate of the central star Thus the IRS 
measurement was well centered on the nebula. The measured 
emission line intensities are given in Table 1, after correcting the 
SH measurements by the factor 2.02 and the SL measurements 
by a factor of 2.35, in the column labeled 'intensity'. The Hfi 
flux found from the infrared hydrogen lines ( especially the lines 
at 12. 37iUm) using the theoretical ratios of iHummer & StorevI 
(Il987h . is 6.15xl0-"erg cm^^ s"', which is about 20% higher 
than the total H/? intensity. This indicates that the LH diaphragm 
measured the entire nebula. The measurement of a higher flux in 
the infrared is within the uncertainties of the various measure- 



ments. Two of the features in the table have been identified as 
belonging to the fullerene molecules tCami et al.ii201Q) . 

4.2. Visual spectrum 

There are only a few visual spectra of Tc 1 . The best is the very 
good spectrum reported b v . Williams et al.l (l2008i) . These authors 
measured the nebula at two positions at either side of the cen- 
tral star, but carefully avoiding the star They used a rectangu- 
lar slit 2"x 4". They correct their intensities for an extinction 
found from the Balmer decrement. These corrected intensities 
are shown in the third column of Table 8 for some of the lines 
of interest to us. Noti ce that the [Nem] U ne at /13869 A was 
too weak to measure. iWilhams et al.l ( l2008h do not report the in- 
tensities of any of the H e i lines; for the intensities of t hese lines 
the spectrum reported bv lKinesburgh & Barlowl d 19941) are used. 
These are shown in column 4 of the t able, where the averag e in- 
tensity is weighted to the spectrum of IWilliams et al.l (l2008h . 



4.3. The ultraviolet spectrum 

There are 14 lUE spectra of Tc 1: Three high resolution short- 
wave spectra, seven low resolution shortwave spectra and four 
low resolution longwave spectra. Only a few lines are strong 
enough for a good identification as a nebular line however. In 
addition the spectrum is of low excitation so that the connection 
between the ultraviolet and visual spectra through the He ii lines 
of /1 1640 A and /14686 A cannot be made. However the nebula 
is small enough s o that almost all of its emission is measured. 
iFeibelmanI (Il983h has reported measuring the [O ii] line at A24-7 1 
A which can be used to connect the ultraviolet spectrum to the 
visual spectrum but we feel that the spectra are too noisy to mea- 
sure this line. We find that only two lines are clearly measurable 
on the low resolution spectra. The high resolution spectra shows 
more lines but because they may be interstellar or stellar we do 
not report them here. Table 10 lists our measurements. The ex- 
tinction correction is made using a value of C=0.33 as found by 
IWilliams et all (I2008h . 
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Table 8. Visual spectrum of Tc 1 . 



Ion 



Intensities ' 



Average 



(A) 




(1) 


(2) 


Intens. 


3726 


[On] 


130 




130 


3729 


[On] 


86 




86 


4340 


Uj 




44.1 


44.1 


4363 


[Oiii] 


0.555 


0.45 


0.55 


4471 


He I 




1.1 


1.1 


4861 


H/3 


100 


100 


100 


5007 


[Oiii] 


124 




124 


5192 


[Arm] 


0.0314 




0.314 


5517 


[CI III] 


0.285 




0.285 


5538 


[CI III] 


0.303 




0.303 


5755 


[Nil] 


1.09 




1.09 


5876 


He I 




9.01 


9.01 


6312 


[Sm] 


0.471 




0.471 


6584 


[Nil] 


95.4 


95.1 


95.4 


6717 


[Sn] 


2.20 




2.20 


6731 


[Sn] 


3.50 




3.50 


7135 


[Arm] 


5.75 


5.29 


5.65 


9069 


[Sm] 


13.0 




13.0 


C(H^) 




0.33 


0.40 





t Ref erences; (1) IWilliams et d] ( l2008h . (2) iKingsburgh & Barlow! 

( fl^4[) . 

(:) indicates uncertain values. 

C is the extinction used by the author 

Table 9. UV Spectrum of Tel . 



A 


Ion 




Intensities 


(A) 


(1) 


(2) (1/^8) 


1906 
1909 

2325 


Cm] 
Cm] 
Cii] 


1.45 
1.16 
3.90 


7.6 14.8 
6.1 11.9 

23.4 45.5 



(l)Measured intensity from in units of 10 '" erg cm - s '. 
(2)Intensity corrected for extinction in units of 10"'^ erg cm"^ s"'. 
I/H/3 is normalized to ^8=100 



4.4. Electron density 

The ions used to determine A^e are listed in the first column of 
Table 4 and in a similar manner as in Table 4. the electron density 
is given. It is about 3000 cm"-'. 



4.5. Electron temperature 

Five ions have lines originating from energy levels far enough 
apart that their ratio is sensitive to the electron temperature. 
These are listed in Table 5. An electron temperature of 9000 
K will be used with an uncertainty of less than 500 K. 



4.6. Ionic and element abundances 

The ionic abundances have been determined using equation 1 
above with an electron temperature of 9 000 K and a density of 
3 000 cm"^. The results are given in Table 10, where the first 
column fists the ion concerned, the second column the line used 
for the abundance determination and the third column gives the 
intensity of the line used relative to Hy6=100. The fourth col- 
umn shows the ionic abundances, and the fifth column gives the 
Ionization Correction Factor (ICF), determined with the help of 



the model described below. In all cases when the ICF is greater 
than 1, the principal ionization stage of that element has been 
observed. 



Table 10. Ionic concentrations and chemical abundances in Tc 1 . 
Wavelength in Angstrom for all values of A above 1000, other- 
wise in fj.m. 



Ion 


A 


lnl./U/3 


A'.on/A'p 


ICF 


A^el./A'p 


He+ 


5875 


9.0 


0.060 


7 


>0.060 


C+ 


2325 


45.5 


1.85(-4) 






C++ 


1909 


26.7 


1.71 (-4) 


1 


3.6(-4) 


N+ 


6584 


95.4 


2.58(-5) 


1.4 


3.6(-5) 


0+ 


3727 


130 


1.87(-4) 






0++ 


5007 


124 


6.8(-5) 


1 


2.6(-4) 


Ne+ 


12.8 


37.5 


6.2(-5) 






Ne++ 


15.5 


1.46 


1.1 (-6) 


1 


6.3(-5) 


S+ 


6731 


3.50 


1.1 (-6) 






S++ 


6312 


0.47 


1.7(-6) 






S++ 


18.7 


14.0 


1.58(-6) 






S+3 


10.5 


0.47 


3.8(-8) 


1 


2.8(-6) 


Ar+ 


6.99 


32.8 


4.2(-6) 






Ar++ 


8.99 


6.5 


7.0(-7) 






Ar++ 


7135 


5.75 


7.0(-7) 


1 


5.1 (-6) 


C1++ 


5538 


0.303 


5.9(-8) 


1.6 


9.4(-8) 


P++ 


17.9 


0.735 


1.67(-7) 


1.2 


2.0(-7) 


Fe++ 


22.9 


0.295 


1.1 (-7) 


1.4 


1.54(-7) 



Intensities given with respect to Hy6=100. 



5. He2-108 



He2-108 (PN G316.1+08.4, also known as IRAS 14147-5156) 
is morphologically quite similar to NGC 1535 and Tc 1. It is 
rough ly circular and has a size at the 10% level of 13.6"x 
12.3" (lTvIenda et al.ll2003h. A somewhat smaUer diameter (11") 
is given bv'A cker et al.l (11992). The size is small enough so that 
most of the radiation can be measured in the lUE diaphragm. 

The 6cm cont inuum radio flux density has bee n measured 
bv lJVrilne&Aiie^ ^975.) as 33 mJv. Itvlilne feAlie^ d 1 98Z) have 
measured the 2cm radio flux density as 43 mJy, which corre- 
sponds to a value of 49 mJy at 6cm. An uncertain average value 
of 39 mJy at 6cm corresponds to an H/3 flux of 1.26 x 10 '^ 
erg cm~^ s~^ The mea sured Hfi flux is 3.7 x 10"'^ erg cm"^ s"' 
(see lAcker et alj[l992l) which leads to an extinction coefficient 
C=0.53. 



5.1. Infrared spectrum 

The IRS measurement of He2-108 was centered at RA(2000) 
14''18'"08.4' and Dec(2000) -52°10'38. 0". This is a slight mis- 
pointing from the center measured by iKerber et al.l (l2003h of 
RA(2000) 14''18'"08.89' and Dec(2000) -52°10'39.7", which is 
presumably the coordinate of the central star. This does not have 
an important effect for the LH measurement because of the large 
LH diaphragm. It does however aff'ect the SH measurement for 
which the Nod 1 measurement measured only part of the nebula. 
The Nod 2 measurement fell within the nebula so that we have 
only used the Nod 2 measurement. By equating the continuum 
measured at 19jum in the LH measurement with the same contin- 
uum measured in the SH Nod 2 we obtain a ratio of LH/SH=2.5. 
The ratio of SH to SL was obtained by equating the [Ne ii] fluxes 
in the two measurements. All fluxes were measured using the 



8 
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Gaussian line-fitting routine. The measured emission line inten- 
sities are given in Table 1, after correcting the SH measurements 
by the factor 2.5 and the SL measurements by a factor of 4.07. 
The Hj6 flux found from the infrared hydrogen lines (especially 
the lines at 7.48yum and 1 2.37;um) using the theoretical ratios of 
dHummer & Storevl[T987l) . is 1.15x10 "erg cm"^ s ', which is 
only slighly smaller than the H/3 found from the radio flux den- 
sity. This indicates that the LH diaphragm measured almost the 
entire nebula. 



5.2. Visual spectrum 

There are two measurements of the visual spectrum of He2- 
108, probably because it is weak a nd not visible to northern ob- 
servato ries. The measurements by 'T orres-Peimbert & Peimberl 
( Il977h are considered by these authors to be less accurate than 
measurements of other PNe they have made. An accuracy of 
about 30% is given by these authors. The visual spectrum has 
also been measured in the 1990 Acker-Stenholm ESO survey of 
southern PNe. The full results of this survey have not yet been 
published, but Acker has sent us the reduced spectrum of this ob- 
ject. The results are given in Table 1 1 . The values are corrected 
for extinction; those in col. 3 by the authors and those in col. 4 by 
us in an attempt to produce the expected Balmer decrement. The 
differences in the blue part of the spectrum reflects uncertainty 
of the measurements. 



Table 11. Visual spectrum of He2-108. 



A 


Ion 


Intensities^ 


Average 


(A) 




(1) 


(2) 


Intens. 


3727 


[On] 


138 




138 


3869 


[Neiii] 


4.7 


20: 


5: 


4340 


Hr 


43 


44 


43 


4363 


[Om] 


1.1 




1.1: 


4471 


He I 


6.8 




6.8 


4861 


HyS 


100 


100 


100 


5007 


[Om] 


200 


167 


190 


5876 


He I 


17 


17.5 


17.5 


6563 


Ho- 


295 


310 


300 


6584 


[Nil] 


105 


90 


100 


7135 


[Arm] 




15 


15 


C(H/?) 




0.4 


0.5 






t References: (1) |T 


orres-Peimbert & PeimbertI (119771). 


(priv.comm.). 










(:) indicates uncertain values. 








C is the extinction (see 


: text) 









(2) Acker 



5.3. The ultraviolet spectrum 

There are 6 lUE spectra of He2-108: 4 low resolution shortwave 
spectra and 2 low resolution longwave spectra. The spectra are 
dominated by the bright central star. Our interest was to look for 
evidence for nebular emission in the carbon ions, either [C m] 
/1 1907 A or [Cn] /12325 A. No [Cii] emission could be seen. A 
[Cm] A19Q1 A line is seen in one of the four shortwave spec- 
tra: SWP14181. Strangely this line is not seen in the other three 
spectra whose exposure time is longer than SWP14181. The ob- 
served line is well above the noise and is at the correct wave- 
length, further the position of all four spectra are the same and 
only the position angle is different. Until new observations are 
available we will regard this measurement as an upper limit to 



the intensity. The observed value is 3 x 10 '''erg cm"^ s ', which 
after correction for extinction becomes 3.3 x 10^'^erg cm"^ s"'. 
There are no other obvious nebular lines in these spectra. 



5.4. Electron density 

There is only a single ion where a line ratio can be used to de- 
termine A^e- This is the [S in] ratio 33.4/18.7 in the infrared. The 
ratio measured, 0.485, leads to a density of 2200 cm"^. The ratio 
of the [Arm] lines in the infrared, 8.99/21.8, is not very depen- 
dent on the density in this range but it is consistent with this 
value of density. The abundances are not dependent on the den- 
sity in this range. 



5.5. Electron temperature 

Only three ions have lines originating from energy levels far 
enough apart that their ratio is sensitive to the electron temper- 
ature. These are listed in Table 5. The temperature found from 
the [O m] lines is uncerta in. A somewhat lower value of 8 730 is 
given by ilVIcKenna et al.l (Il996.) but the authors do not give the 
details of the spectra. We will use a temperature of 8 500 K with 
an uncertainty of about 500 K. 



5.6. Ionic and element abundances 

The ionic abundances have been determined using equation 1 
above with an electron temperature of 8 500 K and a density of 
2 200 cm"^. The results are given in Table 12, where the columns 
are arranged as in Table 12. The ICF is usually 1, except for ni- 
trogen where it is assumed that the ratio N^/N^^= O^/O^^. The 
principal ionization stage has been measured in iron, argon and 
carbon but a small correction has been made for the singly ion- 
ized state. For argon, the similarity of the ionization potentials 
to nitrogen is the basis for the ICF used. 



Table 12. Ionic concentrations and chemical abundances in He2- 
108. Wavelength in Angstrom for all values of A above 1000, 
otherwise in yum. 



Ion 


A 


Int./Bfi 


NUN, 


ICF 


^el./^p 


He+ 


5875 


17 


0.11 




0.11 


C++ 


1909 


<18 


<1.7(-4) 


1.1 


<1.9(-4) 


N+ 


6584 


105 


3.2(-5) 


1.9 


6.0(-5) 


0+ 


3727 


67 


1.4(-4) 






0++ 


5007 


200 


1.35(-5) 


1 


2.8(-4) 


Ne+ 


12.8 


102 


1.72(-4) 






Ne++ 


15.5 


16.5 


1.14(-5) 


1 


2.9(-4) 


S+ 


6731 


2 


1.4(-7) 






S++ 


18.7 


72.2 


7.46(-6) 






S+3 


10.5 


3.6 


4.4(-7) 


1 


8.1(-6) 


Ar++ 


8.99 


25 


2.7(-6) 


1.9 


5.1(-6) 


Fe++ 


22.9 


1.02 


3.6(-7) 


1.2 


4.0(-7) 



Intensities given with respect to 11/3=100. 



6. NGC 6629 



NGC 6629 (PN G009.4-5.0, IRAS 18226-2313) is classified as 
an elliptical, almost round nebula. It is slightly larger than Tc 1 
and He2-108, but smal ler than NGC 1535. Its size down to the 
10% level is given bv JTvlenda et alJ (l2003h as 16.6"x 15.5". It 
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is surrounded by a halo which has a diameter of about 40"and 
is more compressed on the north side. The 6 cm radio flux den- 
sity is given by Griffith et al. ( 1994) as 277 mJy. Milne & Allei 
(11975 ) give a 6 cm flux density of 292 mJy while iMilne & AUei 



Table 13. Visual spectrum of NGC 6629. 



(119821) find a 2 cm flux density of 234 mJy, which corresponds 
to a value of 260 mJy at 6 cm. We will use a value of 270 
mJy for the 6 cm radio flux density which predicts a value of 
HjS=9.3 X lO^^'erg cm"^ s"' for the values of electron temper- 
ature and helium abundance given below. Since the observed 
Hy6=1.18x 10- ' 
orEB-v=0.61. 



erg cm s the extinction constant C=0.896 



6.1. Infrared Spectrum 

The IRS measurement of NGC 6629 was centered at RA(2000) 
18''25"'42.5' and Dec(2000) -23°12'10.1" Thi s is very close 
to the center measured by (iKerber et all l2003h of RA(2000) 
18''25"'42.45' and Dec(2000) -23°12'10.59", which is presum- 
ably the coordinate of the central star. Because the nebula has a 
size close to the sixe of the large LH diaphragm most of the neb- 
ula was within the LH diaphragm. The SH diaphragm measured 
only part of the nebula. By equating the continuum measured 
at 19//m in the LH measurement with the same continuum mea- 
sured in the SH diaphragm we obtain a ratio of LH/SH=2.7. This 
number is somewhat uncertain because the spectrum is rather 
noisy. This ratio can also be obtained by comparing both the 
[Neiii] 15.5/36.0 line ratio and the [Arm] 21.8/8.99 line ratio 
since both of these ratios have only a small dependence in elec- 
tron density and temperature. We obtain an LH/SH ratio of 1.71 
from the [Neiii] lines and a value of 2.44 from the [Arm] lines. 
An average value of LH/SH=2.1 was used. The ratio of SH to 
SL was obtained by equating the [Neii] the [Siv] fluxes in the 
two measurements which leads to SH/SL=1.10. All fluxes were 
measured using the Gaussian line-fitting routine. The measured 
emission line intensities are given in Table 1, after correcting the 
SH measurements by the factor 2.1 and the SL measurements 
by a factor of 2.3. The HyS flux found from the infrared hydrogen 
lines (especially the lines at 7.48iUm and 12 .37/im) using the the- 
oretical ratios of iHummer & StorevI (Il987i) at a temperature of 
8600 K is 8.45xl0-"erg cm~^ s^', which is only slighly smaller 
than the H/3 found from the radio flux density. This indicates that 
the LH diaphragm measured most of the nebula. 



Ion 



Intensities' Average 



(A) 




(1) 


(2) 


Intens. 


3727 


[On] 


36.2 


41.1 


38 


3869 


[Nera] 


40.8 


33.3 


37 


4340 


Hr 


46.3 


46.5 


46.4 


4363 


[Om] 


2.7 


2.8 


2.75 


4471 


He I 


4.6 


4.43 


4.5 


4861 


UJS 


100 


100 


100 


5007 


[Om] 


670 


674 


672 


5517 


[Clm] 


0.45 




0.45 


5538 


[Clm] 


0.4 




04 


5876 


He I 


15.2 


12.0 


14 


6312 


[Sm] 


0.5 




0.5 


6563 


Hff 


286 


289 


287 


6584 


[Nn] 


10.8 


11.0 


10.9 


6717 


[Sn] 


0.6 


0.53 


0.53 


6731 


[Sn] 


0.7 


0.71 


0.71 


7135 


[Arm] 


12.7 


12.1 


12.4 


9532 


[Sm] 


29.2 




29.2 


CiUfi) 




0.8 


0.96 





t References; (D lMilingo elai] ( l2002af) . (2') lAller& KevesI ( 119871) 
(:) indicates uncertain values. 
C is tlie extinction (see text) 



All spectra appear to be dominated by the bright central star. 
This indicates that much of the nebula is within the diaphragm 
but it is difficult to specify exactly how much of the nebula is 
being measured. The only clearly nebular emission is the [C iii] 
/1 1907 A line which is seen in all three shortwave spectra. It 
has the same intensity in all spectra: 2.5 x lO^'-^erg cm"^ s"'. 
Corrected for extinction this becomes 2.2 x lO^^erg cm"^ s"'; 
thus the ratio of the line to H^Sis 23.7 (when Hj8=100). This value 
will be used when determining the carbon abundance but it is a 
lower limit since some of the nebular emission may be outside 
the nebula. Possible [C n] /12325 A cannot be seen. 



6.4. Electron density 

The ions used to determine A^e are listed Table 4. The electron 
density appears to be about 2000 cm"^. 



6.2. Visual spectrum 

There are several measurements of the visual spectrum of 
NGC 6629 in the literature. The two mea surements quo t ed here 
are probably the best. The se are those of iMihneo et al.l (l2002ah 
and Aller & Keves (1987) and are given in Table 13. In addi- 
tion Kingsburgh & English ( 1992) have measured the [O ii] ratio 
3726/3729 to be 1.63+0.50 and the [Sir] ratio 6717/6731 to be 
0.74+0.04. 

The extinction coeficients, C, found from the Balmer decre- 
ment and given at the bottom of the table, are essentially the 
same as that found from the radio flux density (given at the be- 
ginning of this section). 

6.3. The ultraviolet spectrum 

There are seven lUE spectra of NGC 6629: three low resolution 
shortwave spectra and four low resolution longwave spectra. The 
diaphragm was centered near the edge of the nebula for some 
of the spectra and closer to the central star for other spectra. 



6.5. Electron temperature 

Four ions have lines originating from energy levels far enough 
apart that their ratio is sensitive to the electron temperature. 
These are listed in Table 5. An electron temperature of 8 700 K 
is found with an uncertainty of less than 300 K. No temperature 
gradient is apparent. 

6.6. Ionic and element abundances 

The ionic abundances have been determined using equation I 
above with an electron temperature of 8 700 K and a density of 
2 000 cm"^. The results are given in Table 14, where, as in Table 
6, the first column lists the ion concerned, the second column the 
line used for the abundance determination and the third column 
gives the intensity of the line used relative to Hy6=100. The fourth 
column shows the ionic abundances, and the fifth column gives 
the Ionization CoiTection Factor (ICF), determined empirically. 
In all cases but one, when the ICF is greater than 1 , the princi- 
pal ionization stage of that element has been observed. The sin- 
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gle exception is nitrogen, where it is assumed that N^/N^ 

/0++. 



=0^ Table 15. Parameters representing the final model for Tc 1 . 



Table 14. Ionic concentrations and chemical abundances in 
NGC 6629. Wavelength in Angstrom for all values of A above 
1000, otherwise in //m. 



Ion 


A 


Int./HyS 


^ion/Wp 


ICF 


^el./^p 


He+ 


5875 


14.0 


0.094 


7 


0.096 


C++ 


1909 


23.7 


1.9(-4) 


1.1 


2.1 (-4) 


N+ 


6584 


10.9 


3.23(-6) 


14 


4.5(-5) 


0+ 


3727 


38 


3.5(-5) 






0++ 


5007 


672 


4.4(-4) 


1.1 


4.8(-4) 


Ne+ 


12.8 


9.95 


1.8(-5) 






Ne++ 


15.5 


89.1 


6.06(-5) 






Ne++ 


3869 


37 


7.2(-5) 


1 


8.4(-5) 


S+ 


6731 


0.7 


5.3(-8) 






S++ 


18.7 


18.2 


1.77(-6) 






S+3 


10.5 


14.4 


3.9(-7) 


1 


2.2(-6) 


Ar++ 


8.99 


12.7 


1.35(-6) 






Ar++ 


7135 


12.4 


1.67(-7) 


1.2 


2(-6) 


cr+ 


5538 


0.4 


9.7(-8) 


1.2: 


1.2(-8): 


P++ 


17.9 


0.61 


4.5(-8) 


1.2: 


2.0(-7): 



Intensities given with respect to Il/?=100. 



Parameter 


Value 


Ionizing source 




Model atmosphere 


WMbasic 


Teff 


34,700 K 


Logg 


3.30 


Logz 


-0.3 


Luminosity 


1480 Lg 


Nebula 




Density profile 


constant density 2850/cc 


Abundance 


H He C N 




12.000 10.916 8.674 7.590 




Ne Mg Si 




8.431 7.481 5.5 5.778 




PS CI Ar 




5.3 6.203 4.963 6.478 


Size 


6.275" (radius) 


Distance 


1.80 kpc 


Dust grains 


Graphites of single size l.Oyum; 


inner radius 


1.077el6cm 


outer radius 


1.690el7cm 


Filling factor 


1.0 



6.7. Errors 

We refer here to possible abundance errors in all of the PNe 
studied here. This is difficult to specify because there are errors 
due to the measurements, the electron temperature and the ICF. 
There is only a neglible error due to uncertainties in the elec- 
tron density. The measurement error depends on the strength of 
the line; for the stronger lines it is probably less than 15%. The 
values of Tg appear to be independent of the ionization poten- 
tial in all PNe considered here. For at least two of the nebulae 
the uncertainty could be as large as 1000 K. The temperature 
uncertainty plays only a small role for the infrared lines but is 
much more important for the ultraviolet lines. Taken together 
we estimate that for all elements except carbon the abundance 
uncertainties are not more than 20-30% for those elements for 
the ICF is close to unity. When the ICF is higher than 1 .5 the 
abundance uncertaintyies are about 50%. For carbon the ICF is 
usually unity but the abundance is very temperature dependent. 
The error is probably slightly higher, of the order of 50% for this 
element. The largest error for helium occurs in the PNe with low 
temperature central stars where neutral helium is present. This 
occurs in Tel, but may also occur in NGC 6629 and He2-108. 
Other sources of erroe for the helium abundance are probably 
small. 



7. Model 

In order to obtain as nearly a correct model as possible, the star 
as well as the nebula must be considered. Modeling the nebula- 
star complex will allow characterizing not only the central star's 
temperature but other stellar parameters as well (ie., log g and 
luminosity). It can determine distance and other nebular proper- 
ties, especially the composition, including the composition of 
elements that are represented by a single stage of ionization, 
which cannot be determined by the simplified analysis above. 
This method can take the presence of dust and molecules into 
account in the nebular material, when there is any there, making 



it a very comprehensive approach. While the line ratio method 
is simple and fast, the ICFs rest on uncertain physics. To this 
end, modeling serves as an effective means, and the whole set of 
parameters are determined in a unified way, assuring self con- 
sistency. Also, in this way one gets good physical insight into 
the PN, the method and the observations. Thus, modeling is a 
good approach to an end-to-end solution to the problem. We use 
various models in Cloudy's library. 

7.1. Tc 1 

7.1.1. Assumptions 



iTvlenda et all (l2003h give a diameter of 12.9" x 12.2" for this 
PN. We have used a diameter of 12.55" in our modeUng. 

7.1.2. Model results 

Numerous models were run and we found that there was the pri- 
mary problem of fixing the stellar effective temperature. While 
the ionization of carbon, neon and argon indicated a somewhat 
lower Tgff, the observed [O iii] lines required higher value. We 
tried a range of temperatures, distances, densities, density pro- 
files and various model atmospheres. Some models have been 
tried with stellar wind and some without while some with sim- 
ple black body atmospheres for the CSPN. Observed [O iii] line 
fluxes seem to be unusually high and in trying to match them, Ne 
and Ar moved up to higher stages of ionization. Many observed 
lines would suggest a cooler Tgff than what our final model in- 
dicated. In the final model, shown in Table 15, we have used a 
low metallicity model atmosphere mainly to take care of infusing 
more photons in the wavelength region below 912 A without in- 
creasing Tgff. The metallicity is only half that of Sun. The value 
of gravity was also kept as low as possible for a similar boost in 
the input stellar radiation. We note that the galactic latitude of 
the CSPN is only around -9 degrees and the fact that we were 
forced to use a low metallicity model atmosphere for such a low 
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Table U 


». The emission line fluxes (H/3 - 100) for Tc 1 








Label 


Line'' 


Model flux 


Obsd. flux 
(dereddened) 


Label 


line 


Model flux 


Obsd. flux 
(dereddened) 


TOIL 


4861A 


100.00 


100.00 


S3 


6312A 


0.40 


0.46 


C2 


1335A 


10.35 


11.43 


O 1 


6363A 


0.78 


0.05 


C3 


1478A 


0.00 


6.65 


N2 


6548A 


32.43 


30.10 


CI 


1561A 


0.34 


1.26 


N2 


65 84 A 


95.70 


93.20 


C2 


1761A 


0.18 


2.08 


Sll 


6716A 


2.38 


2.15 


C3 


1907 A 


16.56 


14.76 


Sll 


6731A 


3.70 


3.41 


C3 


1910A 


11.85 


11.94 


Ar3 


7135A 


15.52 


6.99 


TOIL 


2326A 


86.80 


44.59 


Oil 


7323A 


9.04 


5.43 


Oil 


3726A 


171.42 


131.10 


Oil 


7332A 


7.23 


4.57 


Oil 


3729A 


91.40 


86.59 


Ar3 


7751 A 


3.74 


1.66 


Sll 


4070A 


0.84 


0.62 


Fe2 


8617A 


0.02 


0.01 


Sll 


4078A 


0.27 


0.18 


C 1 


8727A 


0.05 


0.01 


Fe2 


4244A 


0.00 


0.02 


S3 


9069A 


7.64 


12.51 


Fe2 


4359A 


0.00 


0.01 


TOTL 


9850A 


0.54 


0.04 


TOTL 


4363A 


0.53 


0.56 


Ar2 


6.980m 


14.12 


32.84 


P2 


4669A 


0.02 


0.01 


Ar3 


9.000m 


16.23 


6.49 


03 


4959A 


40.50 


41.71 


S4 


10.51m 


0.51 


0.47 


03 


5007A 


121.92 


123.00 


Ne2 


12.81m 


19.68 


37.41 


Ar3 


5 192 A 


0.08 


0.03 


Ne3 


15.55m 


3.24 


1.46 


Nl 


5 198 A 


0.16 


0.02 


P3 


17.89m 


0.72 


0.73 


Nl 


5200A 


0.06 


0.02 


S3 


18.67m 


12.47 


14.03 


CI 3 


5518A 


0.27 


0.28 


Ar3 


21.83m 


1.05 


0.43 


CI 3 


5538A 


0.29 


0.30 


Fe3 


22.92m 


0.27 


0.29 


Ol 


5577A 


0.04 


0.00 


S3 


33.47m 


4.77 


6.21 


N2 


5755A 


1.23 


1.08 


Si 2 


34.81m 


0.32 


0.30 


Ol 


6300A 


2.45 


0.12 











Absolute HP flux Model: 5.83 x IQ-^'ergs cm-- s'' Obsn: 6.15 x 10""erg.s cm.-^ 5"' 

Notes: 'A' in col. 'Line' signifies Angstrom; 'm' signifies yum. In col. 'Label', we have followed the notation used by Cloudy for atoms and ions; 
this will make identifying a line in Cloudy's huge line list easy. Neutral state is indicated by T and singly ionized state by '2' etc., 'TOTL' 
typically means the sum of all the lines in the doublet/multiplet; or it could mean sum of all processes: recombination, collisional excitation, and 
charge transfer Some elements are represented by usual notation as per Cloudy. 



latitude object shows the extreme to which a modeler is driven, 
when faced with anomalous nebular line emission. 

At this stage we were not sure whether any extra source of 
energy was present, as there was no observational clue, but the 
above facts forced us to look at all different possibilities. It might 
as well be that model atmospheres do not realistically represent 
the stellar radiation. Another fact is that the accuracy of the op- 
tical spectra containing the O in lines is quite high, an eiTor of 
only 8% is quoted by the authors for the spectrophotometry. 
Therefore we surmise that something very interesting is happen- 
ing in the formation of O m lines but are unable to hazard any 
guess. But while we did these adjustments in modeling, lines of 
neon, argon and sulphur gave trouble. The sulphur lines were 
improved to get a better fit by adjusting the DR (dielectronic 
recombination) rates since the DR rates are poorly known. The 
final model output spectra are presented in table 16, where it is 
clear that the model fluxes for most low ionization stages (O ii, 
C n, and Ne ii) are too high. 

To reproduce the observed IR dust continua, we found that 
graphites gave a better fit and used them in our modeling rather 
than silicates. The grains included in the modeling were the 
Cloudy's set called Orion distribution which has a bias towards 
larger grain sizes. The match to the observed IR continua is rea- 
sonable (see Fig. 1). So the final model we present is the best we 
could achieve. Overall we feel our exercise raises more pertinent 
questions than answers as this PN seems to throw lots of chal- 
lenges to our current understanding of nebular physics. We are 
inclined to recommend the abundances as determined by the ICF 
method for this PN. A very pertinent point we want to highlight 
is the fact that this PN shows nebular absorption lines too; see 



'Williams et al.l (l2008h and to the best of our knowledge existing 
photoionization codes are yet to have a provision for computing 
the equivalent widths of such lines so that they can also be com- 
pared with observed values. This would make the criteria for a 
good fit (with observation) more stringent. Presently all mod- 
els published till date used only nebular emission lines alone. 
Incorporatingthe formation of absorption lines would be the next 
major paradigm shift in photoionization modeling. 

7.2. Modeling NGC 1535 

We now describe our unsuccessful efforts to make a photoion- 
ization model for this PN. It is necessary to go into the details 
since it reveals insights into aspects of nebular physics which are 
normally taken for granted as well known. 



7.2.1. Assumptions 

The appearance of NGC 1535 from an image taken by 
ISchwarz et al.l (Il992h is nearly spherical. The IR and radio mea- 
surements described earlier give an absolute H/j flux that is con- 
sistent with a diameter of 45", that includes a low density outer 
zone beyond a high density innerzone of diameter of 19". We 
wanted to include a density profile for the nebula to describe 
the variation of the number density N(H) with the radius and 
so derived a template profile from the //„ image taken from 
http://astro.uni-tuebingen.de/groups/pn/ The image was taken 
on an as-is basis and imported into the IRAF and using cross- 
cuts at different azimuths an average profile was generated and 
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then normalized to have a peak value of 1300/cc. The nebular 
radius was normalized to 22.5", to include the low density re- 
gion as well. We used this profile as a starting value but later 
have experimented with modifications to it as well as tried sim- 
ple constant density models. Though the presence of dust and H2 
molecules is observed in this PN we did not include them in our 
simulations. 



7.2.2. Modeling and its implications 

Though we had tried a series of numerical model computations 
for this object by choosing various options in the parameter 
space, the primary road-block we faced was this: the observed 
fluxes of [Oiv] 25.88/im and the Hen 1640A were too high to 
reproduce while those of specie On, S n and CI m were too low 
to be r eproduced by the models. Many years ago. I Alter & KevesI 
(Il987h modelled this object and he modified the incident en- 
ergy spectrum to get a good fit. Later I Adams & KoppenI (Il985h 
(hereafter AK) had to introduce energy from stellar wind at a T^ 
of 3 X 10^ K to adequately reproduce the observations. While 
Aller's adoptation is purely ad hoc, that of AK is not viable 
physics-wise. AK had predicted that such an extra source of en- 
ergy should make the PN an x-r a y emitter but this is not the 
case as shown by iGuerrero et al.l (12000). Secondly they claim 
that this additional energy source produces a good model match- 
ing observations. But looking at their final results, (see table 
13 of their paper) we find that they did not match the diagnos- 
tic hnes of O n and S n well. They did not include the optical 
lines of [Cliii], Hen 1640A and [Oiv] 25.88jum among others, 
whereas we used a comprehensive set of multi-wavelength spec- 
tral observations in this work. A more important point is that 
they claim that a proper nebular model with rigourous radiative 
transfer would match the observations very well. On the other 
hand, we have tried models having proper windy model atmo- 
spheres (for the CSPN) with Cloudy (wherein radiative transfer 
is handled by escape probability), and it did not work out well. 
We have even tried to infuse additional photons below 226A in 
the input stellar radiation but this did not work out. It is clear 
that AK's suggestion of wind plasma as the extra energy source 
is ruled out. 

As mentioned above Hen and [Oiv] lines are very strong 
and to produce them we needed a model atmosphere with a T^fj 
of 120,000 K. But all other lines then do not match properly 
with such a model. We experimented with a whole range of 
values of the CSPN parameters but did not succeed. We note 
that fro m FUV observations of the continua of the CSPN by 
FUSE, Herald & Bianchi| (|2004|) (HB hereafter) obtain a T.fj 
of only 66,000 K. When we ran a model with the CSPN parme- 
ters as given by HB, it did not reproduce the nebular spectra 
well. Diagnostic lines did not match; He 11 and [O iv] lines were 
weak. More importantly HB give a luminosity of around 4000 
Lq which when used in our model, gave the transmitted flux as 
nearly 80% of the incident flux. We found that a luminosity of 
550 Lq would be sufficient to produce the absolute Hp flux cor- 
rectly. But only a high luminosity of 4000© is compatible with 
the observed FUV continua of the CSPN. This creates another 
problem because this PN is quite complex as it shows hydro- 
gen molecular absorption lines in its FUSE spectrum. These H2 
molecules are attributed to be circum-nebular rather than inter- 
stellar by HB. It would be very difficult to imagine the presence 
of these molecules when the nebula leaks the incident stellar 
radiation enormously in the FUV and UV wavelengths! They 
would simply be destroyed by such a strong radiation. The dust 



[Tc 1 ] Dust cont. emissic 



Wavelength (micron) 



Fig. 1. The IR dust continua of Tc 1. The asterisks represent the 
observation from Spitzer and the continuous curve is the model 
output. 



content is not very high in this PN as shown by HB and so the 
idea of dust grains shielding H2 molecules from the strong ra- 
diation is also ruled out. In our modeling experiments we have 
even tried including a binary CSPN by way of introducing two 
sources of input radiation, a hot and a cooler star. We exper- 
imented with different combinations of temperature and lumi- 
nosity but were not able to reproduce the observations. 

In summary, the best we could intuitively guess was that 
since the complicated ionization strucuture as demanded by the 
observations of nebular emission lines was impossible to re- 
produce, it is probable that the wind streaming into the nebu- 
lar material is directly injecting highly ionized elements, par- 
ticularly O and He. We are led to think of such a scenario un- 
der the given context. We feel that this PN is quite compli- 
cated as far as its physics and chemistry is concerned. It is im- 
portant that this new idea of 'wind streaming' is observation- 
ally explored to establish its possibility, especially as it may 
be applicable other PNe which are difficult to model, such as 
NGC 2392. Shocks have been suggested as the source of energy 
need ed to increse the ioni zation in NGC 2392 and NGC 1535 
(e.g. iPeimbert et alJ Il995|). But if the shocks re sult in a hot 
plasma as assumed by lAdams & KoppenI (Il985h they are in 
disagreement with the observations as discussed above. Thus 
the suggestion of 'ion streaming' is a good possibility for both 
NGC 2392 and NGC 1535. We conclude that the determination 
of abundances is possible, at present, only by the ICF method for 
this PN. 



8. Discussion 

In Table 17 the abundances for the four PNe discussed above are 
summarized in the first four lines, followed by the abundances 
for five other nebulae which have already been determined. As 
me ntioned in the i ntroduction a ll these PNe have been discussed 
by i Mendez et alj ( Il988l 1199 2*) and form a very homogeneous 
group. Not only are these nebulae excited by bright, rather low 
temperature central stars which are rather far from the galactic 
plane, the cen tral stars all hav e spectra indicating that they are 
hydrogen rich. lMendezI (1199 ll) classifies them all as either 0(H) 
or Of(H). 

As may be seen from the table, the abundances of these PNe 
are rather uniform. The oxygen abundance varies by afactor of 
1 .9 but some of this may be caused by the known abundance gra- 
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dient with distance from the galactic plane. To better judge this 
effect, an approximate value of the distance of the PNe from the 
galactic center is given in col. 10 of the table. The nitrogen abun- 
dance, or better still the N/O ratio also has the same low value 
for all the PNe, with the single exception of NGC 2392. This low 
value of N/O is the same as found in the solar atmosphere, which 
is also shown in Table 17. Also notice that the helium abun- 
dance of these PNe is very similar to the solar value. C ombining 
these abundances with the theoretical determinations of lKarakai 
(l2003i) we obtain the following picture. These PNe originate 
from stars of low enough mass so that the second dredge-up or 
hot-bottom burning have not taken place. These processes would 
have increased both the helium abundance and the N/O ratio to 
values higher than observed. For comparison the last two entries 
in Table 17 show the abundances of PNe whose central stars are 
of higher mass and have clearly undergone hot-bottom burning. 

The ratio of carbon to oxygen (C/O) clearly varies for the 
PNe shown in Table 17. The three lowest values are similar to 
the solar value, about 0.5. The highest C/O ratio, that of IC 418 
is almost four times as high. The PNe with the lowest C/O ratio 
probably originate in stars of similar mass to the sun while the 
higher C/O value originate in somewhat higher mass stars which 
have undergone the third dredge-up. Following iKarakasI (l2003h 
the initial mass of the stars with the low carbon abundance is 
between 1 and 1 .SM©, while substantial carbon will be produced 
betweem 1.75 and 2.5M0. According to Weidemann (2000) the 
first group corresponds to a final core mass of 0.55 to 0.57Mq, 
while the group with substantial carbon will have a final mass of 
between 0.59 and 0.63Mo. 

The observed nebular abundances do not permit a more 
quantitative determination of the stellar masses than that given 
above. We can however, determine the stellar masses which are 
predicted by stellar evolution theory to see if they are consis- 
tent with the masses determined from the nebular abundances 
above. To do this we mak e use of the sum mary of stellar evolu- 
tion calculations given bv lBlocken (Il995h in the form of an HR 
diagram where the time of evolution from the AGB is marked 
on each evolutionary track (his Fig. 12). Lines of constant time 
(isochrones) are also shown in this figure. We have fixed the po- 
sition of each of the central stars being discussed on this figure 
by determining the stellar temperature and age. The stellar tem- 
pe rature is taken frorn t he spectra of the star s using the work 
of'Mend ez eTal] d 19881) : FKudritzki et all (l2006h : lPauldrach et alj 
(2004), and is listed in Col.2 of Table 18. In seven of the nine 
cases these temp eratures are the same as are determ ined from the 
nebula (e.g. see iPottasch & Bemard-Salasll2010h . In two cases 
the spectroscopic temperature is lower: NGC2392, where the 
difference is considerable, and NGC1535, where the difference 
is much smaller The reason for this difference is not yet un- 
derstood, but is probably related to the heating of the nebula as 
described above. 

The age is PN is determined from the observed size and 
expansion velocity of the nebul ae. These quanti ties have been 
taken from the values listed bv .Acker et alj dl992l) and are given 
in Cols. 5 and 6 of the tabl e. When two values of velocity are 
listed bv l Acker et al.l (Il992h . the value for the [N ii] is used since 
this line is formed farthest out in the nebula. The velocity is mea- 
sured in the line-of-sight while the measured size is tangential 
but because these nebulae are nearly round it is expected that 
these values may be combined. Having determined the age of 
the nebula its position on the HR diagram is now fixed. The val- 
ues of luminosity and core mass corresponding to this position 
are shown in cols. 9 and 10 of Table 18. While uncertainties in 
the value of size and velocity may be considerable and the age 



determination only reliable to within 50%, the values of lumi- 
nosity and core mass found are much better determined because 
the isochrones are so closely spaced, i.e. the luminosity and core 
mass have a rather small dependance on the age in the low tem- 
perature range of the evolutionary tracks. 

The core mass found in this way is the same as that deduced 
above from the nebular abundances. Even in the prediction that 
the PNe with the high C/O ratios will have somewhat higher 
core masses appears to be fulfilled. We may say that the core 
masses predicted from the PNe abundances are consistent with 
those found from stellar evolution. 

From the luminosity found (in col. 9) and the assumption 
that the central star radiates as a blackbody with the tempera- 
ture given in col.2 and the angular radius found from the stel- 
lar magnitude and extinction listed in cols. 3 and 4, the distance 
can be computed. These distances are listed in the last column 
of Table 18. In 6 of the 9 cases this distance agree s to within 
15% w it h the statistical distance s given either by ICahn et al] 
(Il992h or lStanghellini et aP (l2008h . In the remaining three cases 
(NGC2392, NGC3242 and IC418) the hsted distance is about 
60% higher. 

9. Conclusions 

With the help of Spitzer infrared spectra the abundances in four 
PNe have been determined. These PNe are all excited by rather 
low temperature central stars and have similar morphological 
and kinematic properties: they are all nearly round and are rather 
far from the galactic plane. We are able to show that these nebu- 
lae have rather similar abundances of helium, oxygen, nitrogen, 
carbon and other elements. The resultant abundances are sum- 
marized in the first four lines of Table 17. We then show that 
five other PNe with low temperature central whose abundances 
have been determined using Spitzer infrared spectra and have 
the same or similar morphological and kinematic properties also 
have the same or similar abundances. 

By comparing these abundances with those predicted by nu- 
cleosynthesis models by 'Karakas| (l2003h it is deduced that these 
PNe originate from stars of initi al mass between IM© and about 
2.5M0, which according to iWeidemannI ( 12000^ correspond to 
core masses of between 0.56Mo and 0.63Mo. These values of 
core masses are compared with those determined from stellar 
evolution theory using the observed temperature of the central 
star and the measured age of the nebula. The core masses thus 
found are consistent with each other Two details reinforce this 
consistency. First, the higher core masses from the evolutionary 
theory are foun d in PNe which have high C/O abundance ratios 
as the models of lKarakasI (l2003h predict. Secondly the distances 
found from the stellar evolution are in general values expected 
from statistical distance scales. 

There are a number of uncertainties which must still be con- 
sidered. First, it is not understood why the central star temper- 
ature measured in NGC 2392 is so much lower than that found 
from the nebula. Second, distances found by some researchers 
ar e different than found here. For example, the distances given 
bv lKudritzki et alj (l200d ) for 6 of the PNe listed are 50% higher 
than we have found. This could be due to errors in their de- 
termination of the stellar gravity from uncertain line profiles. 
Furthermore the expansion distances found for two of the nebu- 
lae (NGC 3242 and IC 2448) are at least a factor of 2 lower than 
we have found here. This should be carefully considered in the 
future. 
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Table 17. Elemental abundance of PNe with far-infrared data in addition to optical and UV data. 



PNe 


He/H 


C/H 


N/H 


O/H 


Ne/H 


S/H 


Ar/H 


Cl/H 


R^ 


Ref. 






xlO-^ 


xlO-^ 


xlO-^ 


xlO-* 


xlO-^ 


xlO-* 


xlO-** 


(kpc) 




NGC 1535 


0.091 


1.6 


0.32 


2.7 


0.54 


!.3 


1.1 


6.0 


9.9 


a 


Tel 


>0.060 


3.6 


0.36 


2.6 


0.63 


2.8 


5.1 


9.4 


6.3 


a 


He2-108 


0.11 


<1.9 


0.60 


2.8 


2.9 


8.1 


5.1 




4.1 


a 


NGC 6629 


0.096 


2.1 


0.45 


4.8 


0.84 


2.2 


2.0 


12: 


6.0 


a 


NGC 6826 


0.10 


4.8 


0.58 


3.95 


1.5 


2.6 


1.4 


8.5 


8.0 


b 


1C418 


>0.072 


6.2 


0.95 


3.5 


0.88 


4.4 


1.8 


12 


8.8 


c 


1C2448 


0.094 


2.7 


0.55 


2.5 


0.64 


2.0 


1.2 




8.0 


d 


NGC 2392 


0.080 


3.3 


1.85 


2.9 


0.85 


5.0 


2 2 


13 


8.4 


e 


NGC 3242 


0.092 


1.95 


1.0 


3.8 


0.90 


2.8 


1.7 


7.0 


8.1 


f 


Solar 


(0.085) 


2.7 


0.675 


4.9 


(0.84) 


13. 


(2.5) 


31. 


8.0 


g 


NGC 6302 


0.17 


0.6 


2.9 


2.3 


2.2 


7.8 


6.0 


34. 


6.4 


h 


NGC 6537 


0.149 


1.7 


4.5 


1.8 


1.7 


11. 


4.1 


24. 


6.0 


J 



References: a) Present paper; b) Surendiranath & Pottasch 2008, A&A 483, 519; c) Pottasch et al. 2004, A&A 423, 593; d) Guiles, S. et al. 2007, 
ApJ 669, 1282; e) Pottasch et al. 2008, A&A 481, 393; f) Pottasch & Bernard-Salas 2008, A&A 490, 715; g) Grevesse et al. 2010, ApSS 328, 179; 
h) Pottasch et al. 1999, A&A 347, 975; j) Pottasch et al. 2000, A&A 363, 767 

Table 18. Prediction of central star mass and distance from evolution theory 



PNe 


T.ff 
K 


m^, 


C 


Vel. 
km/s 


Rad. 


t 
10' sec 


R,/Ro 


L./Lo 


Mass 

Mo 


dist. 
kpc 


NGC1535 


66 000 


12.11 


0.08 


20 


10.5 


5.3 


0.543 


5000 


0.59 


2.0 


Tel 


32000 


11.38 


0.36 


12.5 


5.0 


5,4 


2.17 


4500 


0.57 


2.64 


He2-108 


32000 


12.72 


0.53 


12 


5.5 


9.2 


1.96 


3550 


0.56 


3.7 


NGC6629 


46 000 


12.87 


0.88 


6.5 


7.7 


9.5 


0.95 


3800 


0.57 


1.8 


NGC6826 


48 000 


10.68 


0.07 


11 


12.7 


8.4 


0.935 


4080 


0.57 


1.42 


IC418 


36 000 


10.23 


0.33 


12 


6 


4.0 


1.60 


5000 


0.60 


1.25 


IC2448 


65 000 


14.26 


0.27 


13.5 


5 


6.5 


0.54 


4550 


0.58 


4.05 


NGC2392 


43 000 


10.63 


0.22 


53 


22.4 


1.7 


1.55 


7600 


0.62 


1.8 


NGC3242 


75 000 


12.32 


0.12 


27.5 


19 


4.7 


0.43 


5100 


0.59 


1.75 
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